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ABSTRACT 

We present new measurements of the evolution in the Lyman break galaxy (LBG) population 
between z ~ 4 and z ~ 6. By utilizing the extensive multiwavelength datasets available in the 
GOODS fields, we identify 2443 B, 506 V, and 137 i'-band dropout galaxies likely to be at z as 4, 5, & 
6. For the subset of dropouts for which reliable Spitzer IRAC photometry is feasible (roughly 35% of 
the sample), we estimate luminosity- weighted ages and stellar masses. With the goal of understanding 
the duration of typical star formation episodes in galaxies at z > 4, we examine the distribution of 
stellar masses and ages as a function of cosmic time. We find that at a fixed rest-UV luminosity, the 
average stellar masses and ages of galaxies do not increase significantly between z ~ 6 and 4. In order 
to maintain this near equilibrium in the average properties of high redshift LBGs, we argue that there 
must be a steady flux of young, newly- luminous objects at each successive redshift. When considered 
along with the short duty cycles inferred from clustering measurements, these results may suggest 
that galaxies are undergoing star formation episodes lasting only several hundred million years. In 
contrast to the unchanging relationship between the average stellar mass and rest-UV luminosity, we 
find that the number density of massive galaxies increases considerably with time over 4 < z < 6. 
Given this rapid increase of UV luminous massive galaxies, we explore the possibility that a significant 
fraction of massive (10 11 M ) z ~ 2 — 3 distant red galaxies (DRGs) were in part assembled in an LBG 
phase at earlier times. Integrating the growth in the stellar mass function of actively forming LBGs 
over 4 < z < 6 down to z ~ 2, we find that z > 3 LBGs could have contributed significantly to the 
quiescent DRG population, indicating that the intense star-forming systems probed by sub- millimeter 
observations are not the only route toward the assembly of DRGs at z ~ 2. 

Subject headings: galaxies: formation - galaxies: evolution - galaxies: starburst - galaxies: high 
redshift - ultraviolet: galaxies - surveys 



1. INTRODUCTION 

The detailed study of various classes of distant galax- 
ies has enabled great progress in understanding the 
star formation and mass assembly history of normal 
field galaxies (fo r recent reviews see Hopk ins fe Beacoml 
l2?MMusll2!)08l iWilkins et al.ll2008l h Multi-wavelength 
probes have been particularly effective in revealing the 
co-existence of diverse categories of galaxies with red- 
shifts z ~ 2-3. These include the relatively unob- 
scured star-fo rming 'Lyman break ' galax ies (LBGs, e.g., 
ISteidel et ail 1199a IShaplev et all 120051) . the infrared- 
selected massive 'distant red' galaxies (D RGs, e.g., 
iFranx et all 120031: Ivan Dokkum et al.l 12006ft and heav- 
ily obscured sub-mm galaxies which contain both in- 
tensely star-forming and active components (SMGs, e.g. 
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ISmail et alJl!998t IChapman et alJl20"05h . The collective 
study of these populations has revealed that the redshift 
range 1 < z < 3 is a formative one when the bulk of 
the stars in present-day m assive galaxies was produced 
(jHopkins fc Beacoidl200ll . 

Understanding the inter-relationship between these 
various sources is an important goal and intense ef- 
forts are now underwa y to address this issue (e.g., 
Ivan Dokkum et al.ll200l iReddv et al.ll2008fh A relevant 
aspect of this discussion concerns the assembly history of 
objects observed during the redshift interval 4 < z < 6, 
corresponding to a period only 1 Gyr earlier. Such data 
may provide valuable insight into the connection between 
actively star-forming and passive populations as well as 
define the mode of star formation in typical massive 
galaxies. 

Over the last five years, deep multiwavelength sur- 
veys have resulted in t he discovery of large samples of 
LBGs at z ~ 4 - 6 dBouwens et all 120071). Despite 
early controversies ( Bunke r et al.l l2004t IStanwav et all 
120031 iGiavalisco et al.l l2004at iBeckwith et all 120061 ). it 
now seems clear that the star formation density de- 
clines with redshift beyond z ~ 3. Recent evidence 
also suggests th e characteristic lumi nosity is also fading 
(lYoshida et al.ll2006t IBouwens et aTJl2007t iMcLure et al.l 
I2008D . IBouwens et al.l ([2007) attribute this evolutionary 
pattern to the simple hierarchical assembly of galaxies. 
Unfortunately, because of the transient nature of star for- 
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mation probed by the rest-frame UV luminosity function, 
these studies alone provide only an approximate measure 
of the evolutionary processes occurring during 3 < z < 6. 
Key to testing the mode of assembly of galaxies in this 
early period is additional information on the physical 
properties, such as the associated stellar mass and the 
inferred age of the stellar populations. The availability 
of deep IRAC dat a for many of the HST fields enables 
such an approa ch (E g arni et al.l 2005: E vles et all 120051 : 
lYanet al.ll200l . 

In this paper, we aim to improve our understanding 
of galaxy evolution during the first 2 Gyr by systemati- 
cally tracking the evolving stellar content of star-forming 
galaxies in uniformly-selected samples spanning three 
separate redshift intervals between z ~ 4 and 6. By 
adding the additional physical parameters of stellar mass 
and luminosity- weighted ages, we can break degeneracies 
associated with studies that rely only on the UV lumi- 
nosity function. 

An important issue is how an d when the quiescent sub- 
set of DRGs seen at z ~ 2 — 3 (|van Dokkum et aT1l2006t 
iKriek et al.l l2006h assembled their mass. Contemporary 
models of galaxy formation suggest that the number den- 
sity of massive galaxies in creases continuousl y with cos- 
mic time over 2 < z < 5 ([Bower et al.l [20 06), a picture 
supported by the relatively young ages of massive DRGs 
at z ~ 2 (|Kriek et al J 12006). Alternatively, it is poten- 
tially feasible that the bulk of the most massive systems 
at z ~ 2 formed their mass at much earlier times (e.g., 
z > 5) which in turn would imply little evolution in the 
number density of massive galaxies over 3< z < 6. We 
seek to constrain the formation epoch of quiescent DRGs 
by studying the evolving stellar mass function of star- 
forming galaxies over 4 < z < 6. Moreover, we wish 
to understand whether most DRGs were assembled en- 
tirely in intense, dust-enshrouded star formation episodes 
(probed by higher redshift SMGs), or if a significant frac- 
tion were formed in less vigorous star-forming systems 
traced by higher redshift LBGs. By determining the 
number density of massive galaxies in the LBG phase at 
each redshift interval, we hope to estimate the fraction 
of DRGs whose progenitors formed a significant compo- 
nent of their mass in relatively unobscured systems with 
moderate star fo rmation rates (SFR) typic al of the LBG 
population (e.g.. IShaplev et al.ll200ll 120051 ) . 

A further question of significance in characterizing the 
mode of early galaxy assembly is understanding the dura- 
tion of star formation episodes in high redshift galaxies. 
In the context of models in which gas a ccretion dominates 
galaxy growth at high r edshift (e.g.. iBirnboim fc Dekell 
2003; iKeres et al.l [20051: IFinlator et alJ 120071) leading to 
rapid and steady star formation, one might imagine that 
the LBGs seen at z ~ 4 have evolved relatively smoothly 
since their formation at an earlier epoch, creating stars 
at a near constant rate. If their assembly timescales 
are long enough, we would thus expect similarly lumi- 
nous z ~ 5 and 6 LBGs to be less evolved versions than 
their z ~ 4 descendants with lower stellar masses and 
younger ages. Alternatively, it is conceivable that this 
era witnesses a rapid increase in the intensity of star 
formation in individual galaxies, as h as been predicted 
for the early stages of galaxy growth ([Keres et al.1 120051 : 
IFinlator et al.ll2~007h . In this case, galaxies would grow 



along a locus of points (or "main sequence" ) in the M*- 
SFR plane. Thus, if viewed at fixed UV luminosity, the 
observed stellar populations would not vary significantly 
from one redshift to another. Finally, in contrast to these 
"sustained" star formation histories, we may expect star 
formation episodes to occur on much shorter timescales. 
If the past duration of star formation is sufficiently short 
(< 300 Myr), then each dropout sample would be domi- 
nated by newly emerged systems, and we would not nec- 
essarily expect to see significant growth in the M+-SFR 
plane over the redshift range studied. 

Earlier work in this direction has been promis i ng but 
the datasets have been limited. iDrorv et alJ (2005) 
traced the mass assembly of galaxies to high redshift us- 
ing K-band selected samples. However since their study 
did not include IRAC data, re liable stellar masses could 
only be derived to z ~ 4.5. IVerma et al.l ([2007) used 
HST and IRAC data to compare the stellar mass and 
age distribution of a small, robust sampl e of z ~ 5 LBGs 
to those derived by IShaplev et al.l (|2001h at z ~ 3. Most 
recently, iMcLure et alJ (|2008h computed the stellar mass 
function of LBGs at z ~ 5 and 6 by scaling the UV 
luminosity function by the average mass/light ratio of 
the LBG population at high redshift. While these efforts 
have improved our understanding of the stellar content of 
high-redshift galaxies, none has yet systematically traced 
the evolving stellar populations of galaxies over the full 
redshift range, using large and uniformly selected sam- 
ples with no scaling assumptions. 

A plan of the paper follows. In §2, we describe the 
GOODS data used in our analysis. In §3, we discuss 
the color selection used to identify dropouts, and crite- 
ria used to remove contaminants from our sample. We 
close the section by assessing the evolving surface densi- 
ties of the dropouts. In §4, we consider the mid-infrared 
properties of the dropouts and in §5, we discuss the pop- 
ulation synthesis models used to infer stellar masses and 
ages and comment on the uncertainties in the derived 
properties. In §6, we examine the redshift evolution in 
the stellar masses and ages of dropouts at a fixed UV 
luminosity. Using this information, we discuss the im- 
plications for the star formation histories of galaxies at 
high redshift. In §7, we study the evolving stellar mass 
functions of LBGs over 4 < z < 6 and estimate the frac- 
tion of LBGs that evolve into quiescent massive galaxies 
at z ~ 2 — 3. Finally in §8, we compute the stellar mass 
densities of the B, V, and i'-drop samples. 

Throughout the paper, we adopt a A-dominated, 
flat universe with fl\ = 0.7, 17m = 0.3 and Ho = 
701i7okms~ 1 Mpc" 1 . All magnitudes in this paper are 
quoted in the AB system (jOke fe Gunnlll983fh 

2. DATA 

2.1. The GOODS Fields 

We focus our analysis on the data from the Great 
Observatories Origins Deep Survey (GOODS). Detailed 
descrip tions of the datasets a re available in the liter- 
ature ([Giavalisco et al.l l2004b| ). so we only provide a 
brief summary here. The GOODS-S and GOODS-N 
survey areas each cover roughly 160 arcmin 2 and are 
centered on the Ch andra Deep Field South (CDF-S; 
iGiacconi et~aTll2002l) and the H ubble Deep Field North 
(HDF-N; |Wjmams et alJ [199(1 ). Extensive multiwave- 
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length observations have been conducted in each of these 
fields. In this paper, we utilize optical imaging from 
the Advanced Camera for Surveys (ACS) onboard the 
Hubble Space Telescope (HST). Observations with ACS 
were conducted in F435W, F606W, F775W, and F850LP 
(hereafter B 43R, VjiQfh -Zs^n) tow ard GOODS-S and 
GOODS-N (jGiavalisco et all l2004blh The average So- 
limiting magnitudes in the vl GOODS ACS data (cor- 
recting for the flux that lies outside the 0'.'5 diameter pho- 
tometric aperture) are B4 35 =27.8, V606=28.0, ? 775 =27.2, 
and Zg 50 =27.0. We also make use of U-band observations 
of GOODS-N taken with the Kitt Peak National Ob- 
servator y 4-m telescope w ith the MOSAIC prime focus 
camera (jCapak et aHl2004D and of GOODS-S taken with 
the Wide- Field Imager mounte d on the 2.2m MPG/ESO 
telescope ([Arnouts et al.ll200ll ). 

In the near-infrared, we utilize publicly available deep 
J and K-band observations of GOODS-S (PI: C. Ce- 
sarsky) using the ISAAC camera on the Very Large 
Telescope (VLT). The sensitivities vary across the field 
depending on the effective integration time and seeing 
FWHM. Average 5a magnitude limits (corrected for the 
amount of flux that falls outside of the l'.'O diameter aper- 
ture) are J ~ 25.1 and K s ~ 24.6. Toward GOODS- 
N, we m ake use of a new Subaru K-band mosaic of 
the field (IBundv et alj |2009) using the MOIRCS camera 
(jlchikawa et al.ll2006f) . T he specifics of t he ob servations 
are described in detail in IBundv et al.1 (|2009j). The 5a 
limiting magnitude varies between K 8 = 23.9 and 24.1 
(after aperture correction) across the mosaic. 

Deep Spitzer imaging is available toward both GOODS 
fields with the Infrared Array Camera (IRAC) as part 
of the "Super Deep" Legacy program (Dickinson et al. 
in prep, Chary et al. in prep). Details of t he observa- 
tions have been described in detail elsewh ere (|Evles et al.l 
120051 : lYan et all 120051 : IStark et all l2007h so we do not 
discuss them further here. The 5a limiting magnitudes 
of the IRAC imaging are ~ 25.9 at 3.6^m and ~ 25.5 
at 4.5/xm using 2'.'4 diameter apertures and applying an 
aperture correction (see §2.3). 

In a later section (§3.3), we will revisit the sensitiv- 
ity limits of the GOODS datasets, utilizing completeness 
simulations to determine magnitude limits for our survey. 

2.2. Optical Photometry 

ACS photometry was obtained from the GOODS team 
rl.l catalog 8 which ws generated from the vl GOODS 
ACS reduction. The photometric zeropoints adopted in 
the catalog are 25.653, 26.493, 25.641, and 24.843 for 
the i?435-band, i>606-band, i' 775 band, and Zg 50 -band, re- 
spectively. We have corrected for the small amount of 
foreground Galactic extinct ion using the C OBE/DIRBE 
& IRAS/ISSA dust maps of lSchlegel et"~aTl IpOlf l: for the 
GOODS-S field, the color excess is given by E(B - V) 
= 0.008 mag; in GOODS-N, the color excess is given by 
E(B — V) = 0.012 mag. Colors are computed using mag- 
nitudes measured O'.'50-diameter apertures. As detailed 
further in §3.1, total magnitudes are computed using 
a combination of the aperture colors and MAG_AUTO 
SExtractor parameter. 

2.3. Near- and mid-infared Photometry 
8 available from http://archive.stsci.edu/prepds/goods 
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Near-infrared fluxes were computed on each of the 
sources in the GOODS rl.l catalog. We utilized 1"- 
diameter apertures centered on the source positions in 
the A CS images. The seeing varied across the GOODS-S 
and GOODS-N fields as different tiles were taken over 
many nights, so we determined separate aperture correc- 
tions from unresolved sources for each tile. For the J 
and i\g-band ISAAC images the seeing is typically good 
(FWHM = 0'.'4 — 0'.'5), and the aperture corrections are 
« 0.3 — 0.5 mag, determined from bright but unsaturated 
isolated stars measured in 6"-diameter apertures. Like- 
wise, the MOIRCS GOODS-N if s -band seeing was fairly 
constant at FWHM=0'/5; photometry was again com- 
puted in l"-diameter apertures and the aperture correc- 
tions are typically 0.3-0.7 mag. 

For the GOODS IRAC images, magnitudes are mea- 
sured in relatively small apertures («s 1.5 x FWHM which 
corresponds to a diameter of 2'.'4) to maximize the signal- 
to-noise ratio (S/N). We applied aperture corrections 
to compensate for the flux falling outside the aperture: 
these were ss 0.7 mag for the IRAC 3.6 and 4.5 /xm data, 
as determined from bright but unsaturated point sources 
in the images using large apertures. The PSF of IRAC 
leads to frequent blending with nearby sources, so care 
must be taken to ensure that the photometry is robust. 
We discuss our strategy for dealing with this in §4. 

3. SELECTION OF HIGH-REDSHIFT GALAXIES 

The goal of this section is to compile a robust sam- 
ple of dropouts at z ~ 4, 5, and 6. In §3.1, we iden- 
tify B, V, and ^-dropouts using standard selection cri- 
teria (e.g.jGiayalisco et alj|2004at iBeckwith et alj|2006t 
iBouwens et al.ll2007l ). We "excise low-redshift and stellar 
contaminants from these samples using a combination 
of photometric, spectroscopic, and morphological tech- 
niques (§3.2). In §3.3, we compute the completeness lim- 
its of our data and determine the surface density of our 
samples. 

3.1. Dropout Selection 

Galaxies at z ~ 4, 5, and 6 are selected via the presence 
of the Lyman-break as it is redshifted through the B435 , 
V606, an d i'775 bandpasses, respectively. Selection of Ly- 
man bre ak galaxies at these r e dshifts has now become 
routine (IStanwav et alj 120031: iGiavalisco et all l2004at 
Bunker et al.| [2004: Bcc kwith et al.l l2006: Bou wens et al.l 
2007f ) . In order to ensure a consistent comparison of our 
samples to these previous samples, we a dopt color crite- 
ria wh ich are identical to those used in IBeckwith et al. 
(120061) and very similar to those used by IBouwens et al" 
(2007). These criteria have been developed to select 
galaxies in the chosen redshift interval while minimizing 
contamination from red galaxies likely to be at low red- 
shift. We detail our color-selection and S/N limit criteria 
for selecting dropouts below. 

As mentioned in the previous section, galaxies are se- 
lected from the GOODS version rl.l ACS multi-band 
source catalogs. In these catalogs, source detection has 
been performed using the z 850 -band images. 

Conditions for selection as B435-drop : 

B435 - Vera > (1.1 + ^606 - 4so) (!) 
B436 - Veoa > 1.1 ( 2 ) 

^606 -4s <l- 6 ( 3 ) 
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S/N(V 606 )>5 (4) 
S/iV(i 775 )>3 (5) 

In addition we also examine each B-drop for a detection 
in the U-band images described in the previous section. 
B-band dropouts at z > 3.5 should not show strong de- 
tections in the U-band as it lies shortward of the dropout 
filter. We look at each source with a U-band detection 
(2a or more) to make sure that the flux is not spurious or 
associated with another object. In total, we excise four 
sources based on the presence of U-band flux. 

Conditions for selection as V606~drop: 

U 606 - i' 775 > 1.47 + 0.89(z 775 - 4 50 ) 



or 2 (6) 

Veoe - *776 > 1-2 (7) 

<77B-*860<l-3 (8) 

S/N(z' 850 )>5 (9) 
S/N(B 435 )<2 or 

#435 - «775 > V 606 - i' 775 + 1 (10) 

Conditions for selection as i' 775 -drop: 

z 775 - z' 85Q > 1.3 (11) 

S/N(z' S50 )>5 (12) 
S/N(V m& )<2 or 

Vfeoe - z'esa > 2.8 (13) 



In order to utilize these dropout criteria, we must com- 
pute accurate colors and total magnitudes for each source 
in the GOODS catalogs. The colors are computed using 
the aperture magnitudes discussed in §2.2. We adopt the 
i 775 -band MAG_AUTO magnitude as the total z 775 -band 
magnitude for the B435-drops. Total magnitudes are sub- 
sequently computed in the B435, V606j and z 850 -bands 
using the B 435 - i' 775 , V 606 - i 775 , i' 775 - 4so colors mea- 
sured using aperture magnitudes (O'.'50-diameter). Sim- 
ilarly, total magnitudes are computed for the V606 and 
i' 775 -drops taking the z 850 -band MAG_AUTO magnitude 
and B435 - Zg 50 , V 606 - z' 850 , *775 ~ z 850 colors. The S/N 
limits stated above are determined from the measured 
flux and error in the total magnitudes and vary depend- 
ing on the noise background at the location of the source. 
On average, these limits are well represented by the mag- 
nitudes presented in §2.2. 

We examine each source in the dropout catalogs, re- 
moving diffraction spikes from bright stars, image arti- 
facts, and spurious features (primarily along the edges of 
the GOODS mosaic), leaving a sample of 2819 B-drops, 
615 V-drops, and 166 i'-drops. A histogram of the mea- 
sured optical magnitudes for the three dropout samples 
is presented in Figure 1. In the next section, we seek 
to further refine these samples, identifying and removing 
stellar and low-redshift interlopers. 

3.2. Removal of Stellar and Low-z Contaminants 

The Lyman-break selection is well known to in- 
clude a variety of interlopers in addition to the de- 
sired high-redshift galaxies. Standard contaminants in- 
clude cool st ars and dusty or old galaxies at lower- 
redshifts (e.g..|Stanwav et al.ll2004tlBeckwith etalll2006t 
iBouwens et al.ll2006|) . We~first seek to morphologically 
identify and excise stellar interlopers from the dropout 
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Fig. 1. — Histogram of i 775 and z' s50 magnitudes for 2819 B- 
drops, 615 V-drops, and 166 i'-drops in GOODS-S and GOODS-N. 
In each panel, the open histogram corresponds to the total sample 
of dropouts while the filled histogram corresponds to only those 
sources that are isolated in Spitzer (see §4.1). 
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catalogs. To assess the fraction of dropouts at low red- 
shift, we first mine the existing spectroscopic surveys for 
known redshifts of our sample. We then proceed to com- 
pute photometric redshifts for all the remaining dropouts 
without spectroscopic redshifts and compute the uncer- 
tainty in the photometric redshifts. 

We identify stellar contaminants by their morphology, 
removing all bright (z' S5Q < 26.0) sources with SExtrac- 
tor stellarity parameter greater than 0.80. We found that 
this value (very similar to those used in Bouwens et al.) 
was optimal in distinguishing between spectroscopically 
confirmed galaxies and stars. Applying this criteria re- 
moves 51 B-drops (20 in GOODS-S, 31 in GOODS-N), 
23 V-drops (12 in GOODS-S, 11 in GOODS-N), and 8 
i'-drops (3 in GOODS-S, 5 in GOODS-N) from our sam- 
ple. Faintward of this limit, the S/N of the GOODS 
data is too low to reliably identify stars b y their stellarity 
index. Following the approach taken in iBouwens et al.l 
(2006), we estimate the number of faint stellar contam- 
inants that remain in the GOODS dropout catalog by 
examining dropout samples in the Hubble Ultra Deep 
Field (HUDF); the exquisite data quality of the HUDF 
enables point sources to be identified with sufficient S /N 
for sources that would be located at the sensitivity lim- 
its of the GOODS data. Using the same dropout cri- 
teria as in GOODS, we select B, V, and i'-drops in the 
HUDF using pub licly-available zgso-band source catalogs 
(|Coe et al.ll2006l) . We find that the fraction of dropouts 
with stellarity indices above 0.8 over the magnitude inter- 
val 26 < 2850 < 28 is very low for the B-drops (0/171), V- 
drops (0/35), and z'-drops (0/21). T hese results support 
the findings of previous studies (e.g. iBunker et al.1 12004 
iPirzkal et aD2005l : IBouwens et afll2006L 120071 ) that stel- 
lar contaminants are most prevalent at relatively bright 
magnitudes (e.g., z$so < 26). 

Dropouts in the GOODS fields have been observed 
spect r oscop i cally with VLT/FORS2 (Vanzcll a et al.1 
l200l l2005l [20081. f2009h Gemini/GMOS Nod& Shuf- 
fle dStanwav et al.l [2004 . 12001 , and Keck/DEIMOS 
dStanwavet al.ll2004n . The results of the VLT/FORS2 
survey are made public in a large datab ase of 1165 to- 
tal sp ectroscopic redshifts in GOODS-S (jVanzella et al.1 
2008). We search the database (version 3.0) for each 
of the dropouts using an 0.5-arcsecond matching diam- 
eter. Thirty-seven B-drops have spectroscopic redshifts 
ranging between z=3.19 and 4.73 with a median redshift 
of 3.7. We remove the one B-drop with a low-rcdshift 
(z=1.54) identification from our sample. Twenty-four 
V-drops have spectroscopic redshifts at z > 4 with a 
median redshift of 4.81. Finally, 20 «'-drops have spec- 
troscopic redshifts at z > 5, while one i'-drop (which we 
subsequently remove) has a low-z identification (2=1.32). 

Photometric redshifts and their associated probabili- 
ties have been determined for each remaining dropout 
by fitting the observed SEDs against template galaxy 
SEDs of various ages, spectral types, and redshifts us- 
ing the Bayes ian Photometric Redshift (BPZ) software 
(Bcm'tcz 2000). We fit the observe d optical through near - 



IR SED against templates from I Coleman et al 



and the starburst templates from iKinnev et al 



1980) 



1996) 





and a very young (10 Myr) starburst template derived 
from the Chariot & Bruzual 2007 (hereafter CB07) mod- 
els assuming an exponentially decaying star formation 
history with r=100 Myr, no dust, and a metallicity of 0.2 



Photometric Redshift 

Fig. 2. — Top: Comparison of photometric and spectroscopic 
redshifts for dropout s with redshifts measured with FORS2 from 
(Vanzclla ct al. 2008). The center dashed line denotes the unity 
relationship, while the two additional dashed lines correspond to 
Az=z sp ec-Zphot=0.25. Bottom: Distribution of photometric red- 
shifts in B, V, and i'-drop IRAC-clean samples, normalized to the 
total sample size. The median inferred redshift is 3.96, 4.75, and 
6.01 for the B, V, and i'-drop samples, respectively. 

Z Q . The code outputs both a Bayesian and maximum- 
likelihood photometric redshift (and associated probabil- 
ity)- 

In Figure 2a, we present the derived maximum- 
likelihood photometric redshifts for the sample of 
dropouts with known spectroscopic redshifts from 
FORS2. The root mean square (RMS) difference be- 
tween the photometric and spectroscopic estimates is 
0.19, 0.25, and 0.18 for the B-drops, V-drops, and 
i'-drops, respectively, excluding one spectroscopically- 
confirmed source with an anomalously low photometric 
redshift. We find that the RMS uncertainty in our BPZ 
photometric reds hifts is similar to tha t from the GOODS 
MUSIC catalog (|Grazian et all |2006[) for the subset of 
dropouts with spectroscopic redshifts. 

We next attempt to eliminate lower redshift contam- 
inants which made their way through the initial color 
cut, removing those sources with photometric redshifts 
that lie below z p hot < 3, 4, 5 in the B, V, and i'-drops. 
In order to avoid removing viable high-redshift sources, 
we only remove those sources whose summed redshift 
probability distribution suggests that they are very likely 
(> 70%) to lie at low redshift. Sources with near equal 
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maxima at low and high redshift can not be confidently 
included or excluded. However these sources are small in 
number. The inclusion or exclusion of this subset with 
ambiguous photometric redshifts does not affect the ob- 
served optical or infrared magnitude distribution. Fol- 
lowing this appraoch, the BPZ photometric redshifts re- 
quire the removal of 258, 47, and 14 (10, 8, and 9%) 
of the B, V, and i'-drop candidates from the catalogs, 
leaving a sample of 2443, 506, and 137 B, V, and i'- 
dropouts across both GOODS fields. In general, the ex- 
cised sources are either marginally detected in the optical 
or very red in their z' 850 — K colors, as would be expected 
for low-redshift sources which satisfy the dropout crite- 
ria. 

In Figure 2b, we present the distribution of photomet- 
ric redshifts from the remaining high redshift samples 
in bins of Az=0.4. The median and standard deviation 
photometric redshifts of the resulting B, V, and i'-drop 
samples are z phot =3.96±0.29, 4.79±0.25, and 6.01±0.25 
respectively. These distri butions are in rough agreement 
with those computed in (|Bouwens et al.|[2~007f ). 

3.3. Surface Densities and Effective Volume 

The GOODS datasets allow sources to be selected 
down to magnitudes fainter than rriAB = 27; however, 
as is apparent in Figure 1, the completeness begins to 
decline brightward of this limit. Hence, to compute the 
correct density of dropouts, we must account for this 
variation in completeness with apparent magnitude. 

To do this, we define a function p(m, z) that repre- 
sents the probability that a galaxy with apparent mag- 
nitude m and redshift z is recovered in our dataset. 
We compute this function for each of the three dropout 
samples by putting thousands of fake galaxies into the 
GOODS images and recreating a photometric catalog for 
the new image using the selection parameters described 
above. The apparent magnitudes of the fake galaxies 
span i' 775 = 21-28 (B-drops), z' sm = 22-28 (V-drops), 
Zg 50 = 23 — 28 (i'-drops) and their redshifts span the 
interval 3.0< z <4.9 for the B-drops, 4.0< z <5.9 for 
the V-drops, and 5.0< z <7.6 for the i'-drops in steps 
of Az=0.1. The sizes of the fake galaxies are chosen 
to be consistent with distributio n of half-light r a dii de - 
rived for z ~ 4 — 6 galaxies in iBouwens et al.l ([2004) . 
The colors of the fake galaxies depend on the galaxy red- 
shift and spectral energy distribution (SED). In order to 
compute these colors, we use the SED of a CB07 model 
with constant star formation history, an age of 100 Myr, 
and no dust as the intrinsic rest-frame SED of the fake 
galaxies. Allowing for a younger age or a color excess of 
E(B-V)=0.1 in the fake galaxies' SEDs alters the com- 
pleteness by roughly 5%, which would not significantly 
change any of our results. The colors are computed at 
each redshift after applying the appropriate IGM absorp- 
tion (e.g. Madau 1995) to the SED. The probability func- 
tion, p(m, z) is then given by the fraction of fake galaxies 
with apparent magnitude, m, and redshift, z, that satisfy 
the dropout color selection criteria defined in §2.4. 

We next determine the volume observed as a function 
of apparent magnitude fo r each of the dropo ut samples 
following the approach of Ste idel et al.l (|1999h using 



V e s(m) = / dzp(m, z) 
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Fig. 3. — Surface densities of dropout samples. Stars, artifacts, 
and low- 2 galaxies have been removed from the initial sample as de- 
scribed in §3.2. We place the remaining sample in magnitude bins 
spanning 0.5 mags and adjust the measured densities for incom- 
pleteness (§3.3). Our samples are sufficiently complete brightward 
of 27th magnitude (vertical dotted lines) to enable reliable den- 
sity measurements. The resulting surface densities are in excellent 
agreement with previous measurements. 



where p{m, z) is the probability of detecting a galaxy 
at redshift z and apparent magnitude m, and dz ^-j is 
the comoving volume per unit solid angle in a slice dz 
over the redshift interval 3.0< z <4.9 for the B-drops, 
4.0< z <5.9 for the V-drops, and 5.0< z <7.6 for the 
i'-drops. 

The effective volumes are subsequently computed as a 
function of apparent magnitude for the dropout samples 
following equation 14. In general, for the dropout sam- 
ples in a single GOODS field, the effective volumes drop 
from 4-5 x 10 5 Mpc 3 at bright magnitudes (toab=23.5) to 
2 x 10 5 Mpc 3 at fainter magnitudes (raAB=27), reflecting 
~ 50% completeness, on average, at the limit of our sam- 
ples. In Figure 3, we show the derived surface densities of 
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our dropouts, corrected using the magnitude-dependent 
completeness p(ra) derived in equation above. The mea- 
sured values are in excellent agreemen t (generally to 
withi n ler) with previous measurements (|Bouwens et al.l 
l2007h at AB magnitudes brighter than 27. Faintward of 
this limit, the completeness in the GOODS data is too 
low to enable reliable density measurements. 

4. MID-INFRARED PROPERTIES OF DROPOUTS 
4.1. Construction of Spitzer-Isolated Subsample 

Mid-infrared photometry is the crucial component for 
estimating the stellar populations of z > 4 galaxies. 
However the low spatial resolution of Spitzer results in 
frequent blending between nearby sources, making accu- 
rate photometry of individual objects difficult. To en- 
sure that objects in our sample were not contaminated 
by neighboring bright foreground sources, we examined 
each dropout by eye and flagged sources that are hope- 
lessly confu sed with flux from bright, nearby sources. 

Following iStark et al.1 (|2007h . we only apply popula- 
tion synthesis models to those objects which are uncon- 
fused in the IRAC images, which we define as having 
no bright (< 24-25 mag), contaminating sources within 
a ~ 2.5 arcsecond radius around the position of inter- 
est. We find that ~ 35% of the dropouts in our catalogs 
are sufficiently isolated to allow accurate photometry. In 
Figure 1, we show the distribution of i 775 and z' S50 magni- 
tudes of the dropouts in our "Spitzer-isolated" subsample 
compared to the full dropout samples. The magnitude 
distribution of this subsample is very similar to that of 
the parent distribution albeit with a slightly greater pro- 
portion of UV bright sources. We correct for the small 
bias this introduces into the stellar mass functions in §7. 

In Figure 4, we present the relationship between the 
IRAC 3.6/^m magnitude and 1775 or 2:350 magnitude for 
each of our dropout samples, foreshadowing the relation- 
ship between the assembled stellar mass and the current 
star formation rate (uncorrected for extinction). It is 
not necessarily obvious that the ongoing star formation 
in any galaxy should bear any relation to the past star 
formation history, yet immediately apparent is a correla- 
tion between the average optical and mid-infrared flux: 
sources that are brighter in the ACS bandpasses are, on 
average, brighter in IRAC. Also noticeable is that the 
median IRAC flux for a given i' 775 or z 850 flux does not 
change significantly over the redshift range probed by 
the three dropout samples. We discuss the implications 
of these trends in more detail in §6. 

4.2. MIPS Detections 

We have put considerable effort into removing low- 
rcdshift contaminants from our dataset. However, as is 
evident from Figure 4, the dropout samples still con- 
tain red objects, some of which are among the bright- 
est sources detected in IRAC. Considering the fact that 
very massive sources (with bright IRAC fluxes) are likely 
much more common at z ~ 2 than at z > 4, it is clear 
that these sources require more scrutiny before proceed- 
ing. 

We can gain some insight into the likely redshifts of 
this population from 24 /jm imaging with the Multi- 
band Imaging Photometer for Spitzer (MIPS) camera 
(Dickinson et al. in prep.; Chary et al. in prep.). At 






Fig. 4. — Mid-infrared vs. optical flux for the B (top panel), V 
(middle panel), and i'-drops (bottom panel). Individual galaxies 
are shown as small red circles; the large black circles correspond 
to the median 3.6^tm magnitude in bins of 1775 or 2850 magnitude. 
The dashed vertical lines correspond to our adopted completeness 
limit for each dropout sample. Those objects that are brighter 
in the optical are on average brighter in the mid-infrared. The 
relationship between mid-infrared and optical flux is very similar 
for each of the three different redshift ranges considered. 

z ~ 2, the MIPS imaging passband probes the bright 
rest-frame 7.7 /jm feature from polycyclic aromatic hy- 
drocarbons (PAHs). As a result, dusty star-forming 
galaxi es at z ~ 2 are c ommonly detected with MIPS 
(e.g., iReddv et all I2006D . Such PAH features would 
not be detected in sources over the redshift range our 
dropouts sample (4< z <6); hence if any of our sources 
are detected with MIPS, it very well may indicate that 
they lie at z ~ 2. 
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In order to determine what fraction of our catalog is 
detected at 24/im, we visually examine the MIPS data 
of each dropout in the Spitzer isolated sample. While 
the total number of dropouts with MIPS detections is 
small (12/800 B-drops, 3/186 V-drops, and none of the 
i'-drops), it is not negligible. As expected, each of the 
sources detected with MIPS is quite red (z' 850 — TO3.6 > 2) 
in addition to being bright in the IRAC bandpasses. 
These sources thus make up a significant fraction of the 
subset of our dropout samples with bright IRAC fluxes 
(10/25, 2/4 of those with m 3 . 6 <23 for the B and V- 
drops respectively) and are hence sure to strongly affect 
attempts to derive the numder density of massive galax- 
ies. 

While we consider these sources as prime low red- 
shift candidates, it is possible that some of these lie 
at z > 4. Cross-correlating the 24/im-detected sub- 
set with our spectroscopic sample, we find that one of 
the three MlPS-detected V-drops has a spectroscopically 
confirmed redshift of z=4.76. Since there are few strong 
PAH features that fall into the 24 /im filter at this red- 
shift, we propose that the 24 /im emission most likely 
comes from a dusty AGN, a conclusion consistent with 
the point-like morphology in the observed optical-frame. 
Importantly, this establishes that not all 24/im-detected 
dropouts are foreground objects. The MlPS-detected 
subsample thus places an upper limit (1.5% for the B- 
drops and 1.1% for the V-drops) on the number of dusty 
low-2; interlopers remaining in our samples. In subse- 
quent sections, we will derive the evolving stellar popu- 
lations of our dropout sample both with and without the 
24/im detected sources. 

5. DERIVATION OF PHYSICAL PROPERTIES 

We infer stellar masses for the dropout sample by fit- 
ting the latest CB07 stellar population synthesis mod- 
els to the observed SEDs. These models include an im- 
proved treatment of the thermally pulsing asymptotic gi- 
ant braiichJTRAGB) phase of stellar evolution, utilizing 
thelMarigo fc Girardil (|2007f ) evolutionary tracks and the 
iWestera et al.l ( 20021 ) stellar library. In this section, we 
first describe the CB07 models and detail the fitting pro- 
cess that we use to infer the stellar masses and ages of the 
galaxies. We examine the effects of the TP-AGB stars on 
the inferred masses and ages of the d ropouts by compar- 
ing th e CB07 output to that from the lBruzual fc Chariot! 
( 2003) population synthesis models (§5.2). We then com- 
ment on the uncertainties in the inferred properties of ob- 
jects in our samples and stack the infrared photometry 
of the faintest sources to achieve a more reliable estimate 
of the average stellar populations of this subset. 



5.1. CB07 Models 

The CB07 models utilize 220 age steps from 10 5 to 
2 x 10 10 yr, approximately logarithmically spaced. For 
each source, we do not include age steps in excess of 
the age of the universe implied by its photometric or 
spectroscopic redshift. We also exclude unphysically- 
low age steps (e.g., below the dynamical timescale for 
LBGs at high-redshift), which we take to be ~ 20 
Myr adopting measur ements of the mean half-light radii 
( Bouwens et al.l [2004D and typical velocity dispersions 
( Erb et al.ll2006f ) of LBGs at high-redshift. 



The CB07 m odels only inc l ude t emplates with the 
IChabrierl ((20031) and ISakoeterl (fl955T) initial mass func- 
tions (IMF). For consistency with previous efforts, we 
choose to focus on the Salpeter IMF in our analysis. 
The effects of utilizing a Chabrier IMF have been dis- 
cussed in detail e l sewhere, includin g our previous work 
(|Evles et al.1120051 : iStark et al.ll2007Tl . 

We use the high resolution (FWHM 3 A) and 1 A pix- 
els evenly-spaced over the wavelength range of 3300 A to 
9500 A (unevenly spaced outside this range). From the 
full range of metallicities offered by the code, we consid- 
ered both solar and 0.2 Z Q models. From several star 
formation histories available, we utilize the constant star 
formation history and exponentially decaying (r) SFR 
models with e- folding decay timescales r=100 and 300 
Myr. We also utilize a model with a linearly increas- 
ing star formation rate, with a slope set to match the 
observed brightening in the UV LF between z ~ 6 and 
4 (0.74 mag at 1500 A over ~ 620 Myr). The CB07 
spectra are normalized to an SFR of lM yr _1 for the 
continuous star formation model. For the exponential 
decay models, the galaxy mass is normalized to M — > 1 
M Q as t — ► 00. 

We also consider the effects of dust on the in t egrat ed 
SEDs, adopting the reddening model of Calzctti (1997) 9 . 
At each age step, the model SED was folded in with each 
of five different reddening curves corresponding to color 
excesses of E(B-V)=0.0, 0.03, 0.1, 0.3, 1.0. 

For each of the galaxies in our sample, the model SEDs 
were interpolated onto each of the filters' wavelength 
scales after correcting the latter to the rest-frame wave- 
lengths using the redshift appropriate for each object. 
Rather than simultaneously fitting the mass, age, dust, 
and redshift, we fixed the redshift at the best-fit pho- 
tometric redshift specified by BPZ. Where possible, we 
instead utilized spectroscopic redshifts. We discuss the 
effects of holding the redshift fixed later in §5.3. To assess 
the uncertainty in the inferred parameters (e.g., mass, 
age, dust) due to redshift error, we perform the SED fit- 
ting at three redshifts: z — z p hot — 0.25, z p hot, 2phot+0.25. 

To account for Lya forest absorption, we decrement 
the flux shortward o f 1216 A (re st-frame) by a factor 
of exp[-(l + z) 3 46 ] (|Madaul 119951) . We further assume 
the flux shortward of 912 A (rest-frame) is zero due to 
Lyman-limit absorption. 

The resulting template SEDs are folded through each 
of the filter transmission curves to produce optical, near- 
IR, and mid-IR magnitudes for each model. Comparing 
the observed magnitudes and their associated errors 10 to 
the derived model magnitudes, we compute a reduced 
X 2 value for each set of parameters (age, E(B-V), nor- 
malization). For individual sources, non-detections were 
set to the la detection limit and the flux error was set 
at 100% (we also consider the average properties of this 

9 The Calzetti reddening is an empirical law given in terms of 
color excess E(B-V) with the wavelength dependence of the red- 
dening expressed as: 

k(A) = 2.656(-2.156 + 1.509A" 1 - 0.198A- 2 + 0.011A" 3 ) + 
4.88 for 0.12/im < A < Om^m, and 

k(A) = [(1.86-0.48A- 1 )/>-0.1)]A + 1.73 for 0.63/im < A < l.OjUm 
and the flux attenuation is given by: 
Fobs(A) = Fo(A)Kr°' 4E ( B - v >k(A). 

10 We increase the observed magnitude errors to 0.1 when they 
are less than this value to account for calibration uncertainties. 
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faint sample via a stacking analysis in §5.2) 

The fitting routine returned the parameters for each 
model which were the best-fit to the broadband pho- 
tometry (i.e., minimized the reduced % 2 ). In addition 
to selecting those parameters which provide the best-fit, 
we also compute the range of ages, normalizations, and 
E(B — V) values that result in an SED with reduced 
X 2 values within Ax^ edticed =l of the minimum value; we 
take these as the 1-a uncertainties associated with each 
parameter. This normalization was then converted to 
the appropriate units and subsequently used to calculate 
the corresponding best-fit stellar and total mass (and the 
corresponding ltr errors) . The conversion factor between 
normalization and mass for a given age and star forma- 
tion history was obtained from the '4color' files in the 
CB07 population synthesis output. 

5.2. The Effect of TP-AGB stars on Inferred Properties 

at £ > 4 

The TP-AGB stars contribute signi ficantly to t he in - 
tegrated rest-frame near-infrared light. iMarastonl (|2005l ) 
has shown that if the contribution from this phase of evo- 
lution is underestimated for galaxies at z ~ 2, the stellar 
masses and ages inferred from models are typically over- 
estimated. Note that this is not necessarily the case at 
z > 4, since at these rcdshifts the first two IRAC chan- 
nels only probe as far redward as the rest-frame optical, 
where the contribution from TP-AGB stars is less dom- 
inant than the rest-frame near-infrared. In this section, 
we investigate the effect of the TP-AGB stars on the in- 
ferred properties of the B, V, and i'-drops by comparing 
the model output of CB07 and BC03 models. 

In Figure 5a, the stellar masses of the B, V, and i'- 
drop samples inferred from BC03 and CB07 are plot- 
ted against one other for solar metallicity. In comput- 
ing the masses and ages, we fix the redshift at the de- 
rived photometric redshift and the color excess is fixed to 
E(B - V)=0. In spite of the addition of TP-AGB stars, 
the inferred masses from BC03 and CB07 are in very 
good agreement. The median stellar mass of the BC03 
i'-drops is less than 1% greater than that for the CB07 
i'-drops. As expected the effect of TP-AGB is more pro- 
nounced at lower redshifts, but nevertheless we still find 
only a 9% difference between the median of the BC03 
and CB07 B-drop samples. 

At subsolar metallicity, the contribution from TP-AGB 
stars is significant (>20-30% of luminosity for 1.3 Gyr 
insta ntaneous burst) in the rest-frame optical (jBruzuall 
2007); hence, we expect the CB07 models to return lower 
masses than the BC03 models for the 0.2 Z Q templates. 
Indeed, we find that in each dropout sample the varia- 
tions between the two models is more extreme than in the 
solar metallicity case (Figure 5b). In the most discrepant 
case (the B-drops), the median mass inferred from the 
BC03 models is ~ 30% greater than that inferred from 
CB07. 

5.3. Systematic Uncertainties in Derived Properties 

Before examining the results in detail, we explore the 
reliability with which we can determine the stellar popu- 
lations of z > 4 galaxies. There is ample discussion of the 
systematic uncertainties inherent in population synthe- 
sis modeling in the literature (e.g.. iPapovich et al.ll200il : 
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Fig. 5. — Comp arison of the model output from 

Bruzual & Chariot (2003) (BC03) and the updated models pre- 
sented in IBruzuaffl2007r TcB07) which now include the contribution 
of TP-AGB stars. The upper panel shows B-drops (blue circles), 
V-drops (green triangles), and i'-drops (yellow squares) for a solar 
metallicity r=100 Myr exponential decay model. The lower panel 
is identical but with sub-solar (Z=0.2 Zq) metallicity. The dotted 
line in both panels indicates the 1-to-l relationship between the 
BC03 and CB07 mass. The discrepancy between the BC03 and 
CB07 models increases at lower redshifts and lower metallicities. 
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120011: lYan et al.l [20051 : lEvles et al.l [20051 : 
[2001 . We direct the reader to these works 



for detailed discussion and only provide a brief discussion 
below on the typical uncertainties in the inferred masses, 
ages, and extinction for our dropout samples. 

First we must consider whether fixing the redshift at 
a single value biases our determination of the mass and 
its associated uncertainty. To test this, we compute the 
best-fit parameters for the small subset of sources with 
spectroscopic redshifts using two different methods: 1) 
fixing the redshift at the photometric redshift determined 
by BPZ and 2) allowing the redshift to float simultane- 
ously with mass, age, and dust. Comparing these to the 
parameters inferred when the redshift has been fixed at 
its true spectroscopic value gives some indication of the 
error that may be introduced by not allowing the redshift 
to float. We find that the best-fit stellar masses inferred 
by method 1 (fixing the redshift at its BPZ photometric 
rcdshfit) are similar to those determined when the red- 
shift is fixed at its spectroscopic value (the standard de- 
viation in the fractional difference between the two mass 
estimates is 17%). Varying the redshift simultaneously 
with other parameters does not result in better agree- 
ment with the masses inferred from the spectroscopic 
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redshifts (the standard deviation in the fractional differ- 
ence between the two mass estimates is 19%). Moreover, 
we find that the uncertainty in the stellar mass that is 
found by perturbing the redshift by Az=±0.25 from its 
best-fit BPZ photometric redshift and allowing the dust 
and age to vary (42% error) is very similar to that de- 
termined by allowing the redshift to float simultaneously 
with dust and age (48% error). These results suggest that 
our method of using fixed redshifts in our SED fitting 
analysis will not strongly bias the mass determinations 
or decrease the inferred mass error. 

We now move on to quantify the uncertainty in the de- 
rived parameters for sources in our dropout samples. The 
SED fitting procedure is clearly most effective at fitting 
sources that are bright enough to be detected with suffi- 
ciently high S/N (~ 10) in the rest-UV and rest-optical. 
For the bright z' 850 — 24.3 dropout presented in Figure 
6 (typical of bright B-drops in our sample), acceptable 
fits (within &Xreduced = ^ °^ * ne best-fit) are found for 
models with ages between 180 and 321 Myr, color ex- 
cesses between E(B-V)=0.0 and 0.1, and stellar masses 
between 2.0 and 2.8xl0 10 M . If we allow the redshift 
to vary by ±0.25 and consider the range of models that 
produce a reduced y 2 within 1 of minimum reduced \ 2 & t 
the central BPZ photometric redshfit, the inferred stel- 
lar masses vary by less than 5% and ages change by less 
than ~ 30%. Variations in the form of the star formation 
history do change the absolute values quoted above (see 
Shapley et al. 2005 for detailed discussion) , with best fit 
stellar masses increasing by 35% and ages reaching 1.4 
Gyr for a constant star formation model. Nevertheless 
the relative stellar masses and ages within our samples 
should be reasonably well determined for equally bright 
dropouts, unless the star formation history varies consid- 
erably within the sample. 

Extracting robust information from the SEDs of faint 
dropouts is more difficult than for the bright sources dis- 
cussed above, largely because the near and mid-infrared 
datasets are not deep enough to detect very faint objects. 
Sources that are undetected in IRAC have a median mass 
of 2x10 s M Q and age of 80 Myr for the r=100 Myr ex- 
ponential decay model, but their SEDs have acceptable 
masses and ages spanning a factor of 10. In general, we 
find that in order to achieve reasonably reliable masses 
for individual objects, we must limit our sample to those 
dropouts with > 10 9 5 M Q . In Figure 7, we present the 
SED of a typical object near this mass limit (best fit 
stellar mass of 5xl0 9 M Q ); acceptable CB07 fits for the 
r=100 Myr exponentially decay star formation history 
are found using models with stellar masses between 2.4 
and 7xl0 9 M Q , ages between 21 Myr and 360 Myr, and 
color excesses between E(B — V) — and 0.3. While the 
data allow a wide range of ages and dust extinctions, the 
stellar mass can be determined reasonably well. 

To obtain an estimate of the average properties of 
galaxies less massive than this limit, we perform a stack- 
ing analysis on the low S/N infrared data. For each 
of the three dropout samples, we select all sources that 
are undetected and unconfused at 3.6 /im, amounting to 
171 B-drops, 29 V-drops, and 17 z'-drops. For each of 
the dropout samples, the median Zg 50 -band magnitude 
of this faint subset is ~ 26.7. Following the techniques 
discussed in I E vies et al.l (|2007f ). we stack the IRAC data 
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Fig. 6. — SED of a B-drop (typical of bright sources in our 
sample) compared with a range of CB07 model SEDs that fit the 
data. Top: Adopting a model with an exponentially decaying star 
formation history with r=100 Myr, we find the SED is best fit by an 
age of 286 Myr and a stellar mass of 2.3 xlO 10 M with E(B-V)=0 
(green curve) . Allowing the dust content to vary, we find acceptable 
fits for stellar masses ranging between 2.0 and 2.8 XlO 10 M Q and 
ages between 180 and 321 Myr. The model with the lower mass and 
age is shown as the red curve (and has E(B-V)=0.1); the model 
that provides the upper bound to the mass and age is shown in 
purple (and has E(B-V)=0). Bottom: Perturbing the photometric 
redshift by Az = ±0.25, we find that the higher redshift solution 
(red curve) provides a significantly less accurate fit to the data 
(Ax^ educed = 3.06) than that provided by the Az = -0.25 SED 
(green curve) and the SED derived using the best fit photo-z. The 
inferred stellar mass in the lower redshift model (2.3xl0 10 M Q ) 
is within 5% of the best-fit value at the best-fit BPZ photometric 
redshift and the inferred age (203 Myr) is within 30%. 

of each of the undetected dropouts (excluding those with 
< 65 ksec of exposure time) using a median combine with 
no weighting to reduce the effect of contamination by 
neighboring sources. We follow a very similar approach 
in stacking the near-IR data, limiting our sample to the 
objects in GOODS-S given the slightly deeper images in 
that field. 

The stacking procedure described above resulted in 
near and mid-IR detections for the B, V, and i'-drops (see 
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Fig. 7. — SED of a B-drop typical of sources at the adopted 10 
Mq stellar mass limit. Adopting a model with an exponentially 
decaying star formation history with r=100 Myr, we find the SED 
is best fit by an age of 320 Myr and a stellar mass of 4.8 xlO 9 Mq 
with E(B-V)=0 (red curve). Allowing the dust content to vary, 
we find acceptable fits for stellar masses ranging between 2.4 and 
6.7xl0 9 M and ages between 21 and 360 Myr. The model with 
the lower mass and age is shown as the green curve (and has E(B- 
V)=0.3); the model that provides the upper bound to the mass 
and age is shown in purple (and has E(B-V)=0). 



TABLE 1 
Photometry of Stacked Images 



m3.6 M m 



m 4.5/im 



_B-drops (z ~ 4) 



27.03(0.15) 26.51(0.15) 26.79(0.31) 27.29 (0.17) 

V-drops (z ~ 5) 

26.7(0.34) 26.55(0.48) 26.89(0.21) 27.71(0.57) 

i'-drops (z ~ 6)_ 



26.69(0.32) 



26.43(0.42) 



27.23(0.33) 



27.64(0.87) 



Note. — The photometric errors of the stacked images arc provided 
in parenthesis next to the measured magnitudes. 



Table 1 for photometry) . The model fit to the composite 
faint B-drop subsample SED is displayed in Figure 8. For 
the exponential-decay (r = 100 Myr) model, the best-fit 
solar metallicity model has a stellar mass of 1.7x 10 8 M© 
(acceptable fits range between 6.2xl0 7 and 3.2x10 s M Q ) 
and age of 48 Myr. We find similar best fitting models 
for the stacked V-drops (1.5x10 s M Q and 25 Myr) and 
i'-drops (2. 6 x 10 s M and 38 Myr). In all cases, the 
SED is actually slightly better fit by a sub-solar metal- 
licity model (Ay 2 c± 0.4- 0.6) in agreement with the re- 
sults of lEvles et al.l (|2007l ): however, the returned stellar 
masses and ages do not change significantly from those 
obtained using the solar metallicity models (Figure 8). 
Importantly, the inferred masses and ages of the stacked 
SEDs are consistent with what we would have obtained 
by merely taking the average of the masses and ages of 
the individual galaxies included in the stack, suggesting 
that in spite of the large uncertainties, we can still use the 
inferred properties of this faint subset and be confident 
that the average properties are robust. 



Fig. 8. — Fit to composite SED of faint B-drop stack. The best- 
fitting exponential decay t=100 Myr model with solar metallicity 
(red SED) has a stellar mass of 1.7x 10 8 M Q and age of ~ 50 Myr. 
Using a subsolar metallicity, Z=0.2Zq, template (green SED), we 
find a marginally better fit to the data (A% 2 = 0.5) with similar 
inferred masses and ages (2.4xl0 8 Mq and ~ 70 Myr). 

6. NEW INSIGHT INTO THE HISTORY OF HIGH 
REDSHIFT STAR FORMATION 

One of the most basic questions in the study of galaxy 
formation is how galaxies grow over cosmic time. Obser- 
vations have now revealed vast populations of vigorously 
assembling galaxies at z > 4, but it remains unclear how 
star formation proceeds in these objects. On one hand 
it is possible that z > 4 LBGs assemble their mass at a 
constant rate, growing continuously between z ~ 6 and 
4. This "steady growth" scenario may be expected if 
high redshift galaxies are constantly fueled by cold gas 
which is in turn rapidly converted into stars. Alterna- 
tively, galaxies may undergo punctuated star formation 
episodes lasting only several hundred Myr or less, such 
that those LBGs seen at z ~ 5, for example, are not 
necessarily related to those seen at later times (e.g., at 
z ~ 4) or at earlier times (e.g., z — 6). This "indepen- 
dent generations" scenario could arise if relatively rare 
events (e.g., major mergers, instabilities) are predomi- 
nantly responsible for sparking luminous star formation 
at z > 4, or alternatively if feedback effects act to shut 
off star formation after several hundred Myr. 

Testing the validity of these scenarios has proven dif- 
ficult with rest-UV observations alone, as these studies 
only probe recent star formation activity, revealing little 
about the physical properties of the galaxies being as- 
sembled. But it is the evolution of such properties (i.e., 
stellar mass, age, dust) which can potentially illuminate 
the relationship between galaxies at z ~ 4 and z ~ 6. 
If the LBGs seen at z ~ 4 have been forming stars at a 
constant rate since z ~ 5 — 6, then the higher redshift 
objects of a given luminosity should appear as "scaled- 
down" versions of those at z ~ 4, with lower masses, 
ages, and reduced dust content. In contrast, if the LBG 
population at z ~ 4 is dominated by galaxies which were 
either quiescent or forming stars less vigorously at z ~ 5 
and 6, then we would not necessarily expect the stel- 
lar populations to show signs of growth between the two 
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different epochs. 

In this section, we attempt to distinguish between 
these various scenarios by exploring the evolution of 
the relationship between the stellar mass and UV lumi- 
nosity (hereafter the M*-Mi5oo relation). For a given 
galaxy, this diagnostic illuminates the ratio of the past 
and present star formation. At lower redshifts, this re- 
lationship has been use d effectively to constrain stellar 
mass growth in galaxies (|Noeske et al.lf2007t iDaddi et al.l 
l2007f) . Studying the redshift evolution of this relation for 
a large sample of high redshift dropouts reveals whether 
the typical specific star formation rate (SFR/M*) and 
stellar mass of UV luminous galaxies changes signifi- 
cantly between z ~ 4 and 6, allowing new constraints 
to be placed on the past star formation history. 

In Figure 9, we present the inferred stellar masses 
of the B, V, and i'-drops as a function of absolute 
magnitude at rest-frame 1500 A (Misoo)- We com- 
pute the absolute magnitude from the observed «'-band 
magnitudes for the B-drops, and the z'-band magni- 
tudes for the V- and i'-drops using photometric redshifts 
where spectroscopic confirmation is not available. Ide- 
ally we would correct the absolute magnitudes for dust 
extinction, but doing so can introduce considerable er- 
ror int o the relation. For e xample, using the Meurer re- 
lation (jMeurer et al.l fl999l ) we find that the extinction 
at 1600 A, Aigoo, is related to the rest-UV color for 
dropouts at z — 4 by the following relation: Ai6oo=0-415 
+ 10.9(i' 775 -Z85o)- This indicates that photometric uncer- 
tainties as small as 0.1 mags in the rest-UV color trans- 
late to more than 1 mag of uncertainty in the extinction 
correction. Further error arises from uncertainty in the 
redshifts and, to a lesser extent, the ages of the dropouts, 
both of which alter the relation between extinction and 
color derived above. Together, these sources of error 
would undoubtedly add significant random scatter to the 
M,t-Mi5oo relation, quite possibly washing out any trends 
that may exist with luminosity or redshift. We thus limit 
our analysis to the relationship between stellar mass and 
emerging luminosity (i.e., the luminosity inferred from 
the flux that escapes the galaxy) and comment below on 
the specific effects that not correcting for dust extinction 
may have on our results. 

For each of the dropout samples, we find that the me- 
dian stellar mass increases monotonically with increas- 
ing emerging UV luminosity (forming a so-called "main- 
sequence"), albeit with very significant scatter in each 
magnitude bin. Si milar relations are commonly seen at 
lower redshifts fe.g-. lDaddi et al.ll2007f) and are predicted 
to ar ise at high redshift in hydrodynamic simulations 
(e.g., iFinlator et al.ll2~006h . In the idealized scenario in 
which the bright dropout population emerges at z ~ 6.5 
and steadily grows in stellar mass via a constant star for- 
mation rate (with no new galaxies emerging as luminous 
dropouts), we would expect the typical dropout stellar 
masses to increase by a factor of ~ 7 in the time between 
z ~ 6 (when the galaxy has an age of ~ 90 Myr) and 
z ~ 4 (when the galaxy has an age of 700 Myr). 

Examining the top, middle, and bottom panels of Fig- 
ure 9, we do not see such a systematic increase in the 
normalization of the M*-Mi5oo relation between z ~ 6 
and 4 (Table 2). While the lowest luminosity bin above 
our completeness limit may show some signs of marginal 
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Fig. 9. — Stellar mass vs. absolute magnitude at 1500A (un- 
corrected for dust extinction) over z ~ 4 — 6. Small red circles 
correspond to inferred stellar masses and rest-UV absolute magni- 
tudes for individual B-drops (top), V-drops (middle) and i'-drops 
(bottom) assuming a r = 100 Myr exponential decay model. The 
dark solid circles are the median stellar mass in each magnitude 
bin. The relationship at z ~ 4 is overlaid on the z ~ 5 and z ~ 6 
panels as a black solid dashed line. The vertical solid lines repre- 
sent the adopted completeness limits for each sample. The median 
stellar mass increases monotonically with M1500 in each dropout 
sample; however at a fixed M1500, the median stellar mass does not 
decrease significantly with increasing redshift, as may be expected 
in simple steady growth models. 
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TABLE 2 

Evolution of Galaxy Properties over 4<z<6 



M uv 


Num 


M* (r=100) 


M* 


(CSF) 


Age 


(mag) 




(10 s M Q ) 


(10 1 


* Ms) 


(Myr) 



-22.5 
-21.5 
-20.5 
-19.5 



15 
73 
227 
370 



-22.5 
-21.5 
-20.5 
-19.5 



3 
20 
(57 



B-drops (z ~ 

144 (31-410) 
77 (24-380) 
29 (4.0-130) 
5.5 (2.0-42) 

V-drops (z ~ 

304 (180-830) 
109 (5.4-234) 
27 (2.2-158) 



4) 



223 
114 
38 
7.3 



180 
203 
181 
143 



5). 



520 
112 
3(5 



286 
227 
181 



-22.5 
-21.5 
-20.5 
-19.5 




8 
11 



_ i'-drops (z ~ 6) 

160 (24-530) 
13 (1.3-90) 



226 
27 



286 
102 



Note. — The stellar masses and ages are inferred from models us- 
ing a Salpctcr IMF and solar metallicity. In column 3, we present the 
median stellar masses (and the range of masses spanned by the middle 
80% of the distribution) determined from an exponentially-declining 
star formation history with r — 100 Myr. In column 4, we provide the 
median stellar masses inferred assuming a constant star formation his- 
tory. In column 5, we present the median ages inferred for the r— 100 
Myr models. 

stellar mass growth between z ~ 6 and 5 (i.e., a x 1.3— 2.0 
increase in the median stellar mass) , this same bin shows 
no growth between z ~ 5 and 4 in spite of the improved 
mass estimates and statistics available at the lower red- 
shifts. The most luminous bins actually show a slight de- 
crease in the average stellar mass with cosmic time. As 
shown in Table 2, the absence of a systematic increase in 
the average stellar masses is not strongly dependent on 
the chosen star formation history. Overall these results 
seem to imply that the ratio of median stellar mass to 
emerging UV luminosity does not evolve significantly for 
LBGs over z = 4 — 6. A galaxy with a given M1500 at 
z ~ 6 will, on average, have the same assembled mass (to 
within a factor of ~ 2) as a galaxy seen at z ~ 4 with the 
same M1500. This suggests that the specific star forma- 
tion rate evolves weakly over 4 < z < 6, indicating that 
the typical duration of past star formation for a galaxy 
of a given luminosity does not vary significantly between 
z ~ 6 and 4. 

While the inclusion of a dust correction would shift 
the M*-Mi5oo relation brightward (i.e., to the left in 
Figure 9), it would likely not lead to an increase in 
the normalization of the M^-Miso n relation over time . 
As has been shown e lsewhere (e.g.. IStanwav et aT]|2005t 
iBouwens et al.l [20071) galaxies at Mi 50 o < —19.8 do po- 
tentially become marginally dustier between z ~ 6 and 4 
which would cause the z ~ 4 M*-Mi5oo relation to shift 
slightly more than the z ~ 5 and 6 relations. Since a 
relative shift brightward is roughly the same as a shift 
toward lower stellar masses at fixed M1500, this would 
actually have the effect of slightly decreasing the nor- 
malization of the M*-Mi5oo relation over the 4 < z < 6 
rcdshift range. 

6.1. Testing the Steady Growth Scenario 

We now attempt to place the M*-Mi5oo relation pre- 
sented in Figure 9 in the context of the steady growth 
scenario discussed at the outset of this section. To ap- 



proximate this scenario, we assume that galaxies follow 
a constant star formation rate. Using solar metallicity 
templates and allowing the dust content to freely vary, 
we find that the typical implied star formation lifetimes 
of the B-drops are in excess of 700 Myr in the two bright- 
est bins (Figure 10), implying that the precursors of the 
majority of B-drops would have been equally luminous 
in V and z'-drop samples. We therefore would have ex- 
pected to see a steady decline in the typical stellar masses 
and ages of UV luminous galaxies at z ~ 5 and 6. The 
fact that the z ~ 6 LBGs are, on average, quite sim- 
ilar in mass and age to those of similar luminosity at 
z ~ 4 suggests that the precursors of most of the z ~ 4 
LBGs were not at the same UV luminosity (or perhaps 
were even quiescent) at z ~ 6. Similar arguments are 
arrived at by considering the observed UV luminosity 
func tion (UV LF) of the d ropouts over this redshift inter- 
val (|Bouwens et al.|[2007l ). The large implied ages men- 
tioned above combined with the constant past star for- 
mation rate suggest that the number density of galaxies 
at a given M1500 should not fall off significantly with in- 
creasing redshift, in sharp contrast to the observed trend. 
Indeed, the distribution of ages implied by this constant 
star formation rate scenario is such that the number den- 
sity of galaxies at Mi 50 o < —21.5 should not decrease by 
more than a factor of ~ 2 over 4 < z < 6, in marked 
contrast to the observed xlO decline in the bright end 
of the luminosity function at M1500 < —21.5 over this 
redshift interval (|Bouwens 

We therefore argue that in order to explain the con- 
stancy in the VL-M1500 relationship and the decline in 
the UV luminosity function, we require a star forma- 
tion model that implies that a significant fraction of UV- 
luminous galaxies at z > 4 emerged at their present 
UV luminosity recently enough such that they would not 
have been visible as dropouts at the same UV luminosity 
just ~ 300 Myr earlier (e.g., the time corresponding to 
Az ~ 1 at 4 < z < 6). With a constantly increasing 
flux of newly luminous galaxies, the galaxy population 
at each epoch would likely be dominated by relatively 
young (< 300 Myr) objects. Hence the main sequence 
of the M*-Mi5oo relation would not necessarily grow sig- 
nificantly with cosmic time as one would expect in the 
steady growth scenario. 

We note that this requirement of a young population 
is not inherently inconsistent with constant, sustained 
star formation in galaxies. Indeed, if each of the z > 4 
dropout populations is dominated in number by newly 
emerged luminous galaxies, it is possible that systems 
which are UV luminous at z ~ 6 remain luminous to 
z ~ 4, but make up such a small fraction (by number) of 
the z ~ 4 LBG population so as not to affect the main se- 
quence of the M*-Mi5oo relation. The key problem with 
this picture, as hinted by our results above, is that con- 
stant star formation models struggle to produce young 
enough ages to be consistent with the need to have each 
epoch (i.e. z ~ 4) dominated by systems that weren't 
present at the previous epoch (i.e. z ~ 5). 

It is valuable to consider the systematic uncertainties 
that could arise following degeneracy in derving both age 
and dust extinction measures from our photometric data. 
If our dust estimates are systematically too low, then the 
ages quoted above are likely overestimates. In addition, 
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Fig. 10. — Ages inferred for B-dropout sample from constant 
star formation scenario as a function of absolute magnitude at 
1500 A (uncorrected for dust). The solid red circles give the me- 
dian inferred age for CB07 models with solar metallicity while the 
blue circles give the median ages for models with 20% solar metal- 
licity. The error bars on each datapoint correspond to the 25th 
(lower bar) and 75th (upper bar) percentile in the distribution of 
ages in each absolute magnitude bin. The horizontal dashed line 
denotes the amount of time spanned between z ~ 4 and 5. Thus 
if the inferred ages are above the dashed line, it indicates a large 
fraction of the B-drops began forming stars beyond z ~ 5, im- 
plying an abundant precursor population with scaled-down masses 
relative to those measured at z ~ 4. The lack of evolution in the 
normalization of the M*-Mi5oo relationship over z ~ 4 — 6 (Figure 
9) casts doubt on this scenario. 

if we adopt lower metallicity templates in our models 
(e.g., 0.2 Z Q ) as may be appropriate for z > 4, we find 
ages that are a factor of ~ 2 lower than those in Ta- 
ble 2 (Figure 10). Thus while the relatively large ages 
implied by the constant star formation rate scenario do 
cause some tension with the observed evolution in the 
M*-M 150 o relation and UV LF, there are feasible steady 
growth scenarios that are potentially consistent with the 
observations. Improved spectroscopic completeness and 
infrared photometry will soon enable us to lessen the un- 
certainty in the dust content and to provide estimates of 
the metallicities, both of which will help test the feasi- 
bility of these scenarios. 

A possible variant on the "steady growth" scenario is 
a situation where star formation grows increasingly vig- 
orous as galaxies grow in mass at early times. While 
such star formation histories are not typically used in 
fitting high redshift galaxies, the y are predicted in SPH 
simulations (|Finlator et al.l l2006f) and have successfully 
been used to infer the properties of two z > 6 galaxies 
(|Finlator e t al. 200^). In this case, galaxies would grow 
along the "main sequence" of the M*-Mi5oo relation of 
Figure 9, which results in little evolution of the normal- 
ization of the M*-Mi5oo relation over 4< z < 6, similar to 
the observed trend. This picture is also able to explain 
the observed luminosity function evolution. If galaxies 
in the luminosity range probed by our study steadily 
brighten by ~ 0.7 mag at 1500 A in the ~ 620 Myr 
between z ~ 6 and 4 (matching the evolution seen in the 
characteristic UV luminosity, MJj V ), then this scenario 
would be fully consistent with the observed evolution in 
the UV luminosity function. 

Importantly, a significant obstacle for each of the 



"steady growth" scenarios described above comes from 
observations of the clustering of dropouts at z > 4 
(|Ouchi et al.ll2004tlLee et al JI2006L l2008f > . By linking the 
dropouts of a given UV luminosity to a dark matter halo 
mass, constraints can be placed on the fraction of halos 
which are UV-bright, from which a star form ation duty 
cycle can be determined. Using this approach, iLee et all 
( 2008) concluded that the typical duration of star forma- 
tion activity in the dropouts at z ~ 4 must last no longer 
than ~ 400 Myr. Hence according to these results, those 
galaxies that are active at z ~ 6 would no longer be ac- 
tive at z ~ 4, in direct conflict with the sustained activity 
inherent to this picture. 

In summary, we find that it is difficult to reconcile 
the observed data (UV LF, stellar populations, and clus- 
tering) with a scenario in which the galaxies observed 
at z ~ 4 have been steadily growing at a fixed UV lu- 
minosity since z ~ 6. Instead the data require z > 4 
dropout galaxies to have been less luminous several hun- 
dred million years prior to the redshift at which they are 
observed. Hence we find that a promising means of rec- 
onciling the steady growth model with the observations 
could be a scenario whereby the star formation rates of 
galaxies are, on average, gr owing more vig o rous a s their 
mass grows, as proposed bv lFinlator et al.1 (|2006l ). 

6.2. Testing the Independent Generations Model 

We now consider the feasibility of the "independent 
generations" model mentioned at the outset of this sec- 
tion. In this picture, the galaxies we observe at z > 4 
undergo shorter timescale star formation episodes lasting 
no more than ~ 300-400 Myr. In order to achieve these 
relatively punctuated episodes of luminous activity, each 
star formation event must eventually decline. We model 
this behavior using the simple exponential decay models 
discussed in §5.3. The typical inferred ages for B-drops 
using solar metallicity models with decay constants of 
r=100 (300) Myr range from 140 (230) to 200 (450) Myr 
at M1500 < —19.5, significantly lower than the constant 
star formation models discussed above. These ages are 
comparable to the time between z ~ 4 and 5, indicating 
that with these assumptions the z ~ 4 dropout popu- 
lation is dominated by sources which were not actively 
forming stars at z ~ 5, as required based on our dis- 
cussion above. Unless the cold gas reservoirs present in 
newly emerged galaxies at a given UV luminosity vary 
strongly between z ~ 4 and 6, we would expect to see 
very little evolution in M^-Misoo relation, similar to the 
observations presented in Figure 9. 

Since the inferred ages of the dropouts in this episodic 
star formation scenario are typically small enough such 
that they would not have been active in higher red- 
shift samples, there is no inconsistency with the observed 
decrease in the number density of luminous sources at 
z > 4. The only requirement is that whatever process 
initiates the star formation activity must steadily pro- 
duce more vigorous episodes over cosmic time. Given 
the expected growth in the underlying dark matter halos 
hosting the dropouts, this is not an unreasonable require- 
ment. 

Unlike the "steady growth" scenario discussed previ- 
ously, the exponentially declining models predict that 
the dropouts will grow progressively fainter over time. 
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Assuming models with decay constants of 100 Myr, the 
sources will become nearly 3 magnitudes fainter in M1500 
in just 300 Myr, and thus would almost certainly escape 
selection in future dropout samples with similar mag- 
nitude limits unless new star formation episodes have 
been initiated. The short (< 500 Myr) star formation 
timescales implied by this scenario are clearly in much 
better agreement with those i nferred from the clustering 
analysis previously discussed (|Lee et al.ll2008f h 

A potential drawback of the exponentially decaying 
models discussed above is that in order to reproduce 
the observed SEDs of the dropout populations in the 
required short timescales, we require the past UV lu- 
minosities to exceed those measured presently. For the 
r=100 Myr models, we would expect the typical galaxy 
to have been a factor of ~ 3 brighter at 1500 A just 
150 Myr earlier (with this very luminous phase lasting 
for a duration of ~ 50 Myr). This implies an abundant, 
precursor population of very young, luminous objects, 
which are seemingly not present in sufficient number in 
the observed dropout sample. This problem could po- 
tentially be explained if this lumin ous phase is generally 
heavily extincted by dust. Indeed, IShaplev et alJ (|2001l ) 
found that the most luminous z ~ 3 LBGs tend to be 
very young and dusty. Further spectroscopy and higher 
S /N near-IR data will soon help test this possibility. 

One prediction from this episodic star formation model 
is that there should be a significant repository of stellar 
mass in objects which underwent a UV luminous phase 
at an earlier epoch but are currently not actively forming 
stars. A careful comparison of rest-UV and rest-optical 
selected samples at z ~ 4 could thus yield significant 
insight into the mode of star formation that is dominant 
at z > 4. By examining the range of star formation 
rates present in galaxies of fixed stellar mass, we should 
be able to distinguish between the steady growth and 
independent generations scenarios. The reliability of this 
test relies on having samples with secure redshifts and 
accurate extinction corrections, requiring both additional 
spectroscopy and improved IR photometry. 

Before concluding this section, we must comment on 
several caveats that apply to our discussion. The first is 
the possibility of an evol ving initial mass function (| Pavel 
2008; van Dokku m 2008) in which the characteristic stel- 
lar mass in creases towa rd higher redshifts. As explained 
in detail in Dave (2008), this would alter the current star 
formation rates and the inferred stellar masses in such a 
way to decrease the normalization of the M*-Mi5oo rela- 
tion at progressively earlier times. Clearly this does not 
help reconcile the observations with the steady growth 
picture. The second caveat is our neglect of the contri- 
bution of galaxy mergers to stellar mass growth. While 
minor mergers are expected to be very common at high 
redshift, they are likely to be gas-rich interactions. Such 
merger events would help provide galaxies with the gas 
required for star formation; however, the established stel- 
lar mass gained in a gas-rich merger would likely not be 
significant compared to the stellar mass growth that is 
occuring due to star formation. This view is supported 
by simulations which predict that stellar mass growth is 
dominated by star formation and not mergers a t the high 
redshifts we are considering (|Keres et al.ll2005l ). 

In summary, the episodic star formation model (ap- 



proximated by a single component exponentially declin- 
ing star formation history with luminous lifetimes of 300- 
500 Myr) can explain the evolution in the M*-Mi5oo rela- 
tion, the UV LF, and the clustering properties at z > 4. 
If this mode of star formation is dominant at z > 4, 
it would signal a potential shift from the mode thought 
dominant at z ~ 0.5 — 2, where the lack of significant scat- 
ter in the dust-corrected SFRs of galaxies at fix ed stel- 
lar mass (|Noeske et al.l l2007t iDaddi et al.l l2007t ) argues 
against the feasibility of sta r formation being dominated 
by short-timescale episodes ()Davell2008l ). However, given 
the minimal spectroscopic completeness and low S/N of 
the IR data at z > 4, it is too premature to rule out 
the sustained mode of star formation for z > 4 galax- 
ies. With new IR imaging and extensive spectroscopic 
surveys for high redshift dropouts currently underway, it 
will soon be possible to improve the fidelity of the re- 
sults presented here as well as to robustly test several of 
the additional predictions of the episodic star formation 
mode described above, both of which should lead to fur- 
ther progress in our understanding of the history of star 
formation at z > 4. 

7. CONNECTING THE Z > 4 MASS ASSEMBLY HISTORY 
TO QUIESCENT GALAXIES AT Z ~ 2 

The assembly history of massive galaxies at high- 
redshift [z > 2) has been a very active area of research 
in recent years. Near-infrared photometry of rest-UV se- 
lected objects at z ~ 2—3 has revealed a si gnificant popu- 
lation of ma ssive galaxies at high redshift (|Shaplev et al.l 
120011 [2005] ). In addition, near-infrared color selected 
samples have proven very effective at is olating large sam- 
ples of DRGs (e.g., iFranxet al.ll2003f) . The red colors 
of these objects likely arise from evolved stellar popu- 
lations or dust (or some combinatio n of the two) . The 
addition of Spitzer photometry (e.g . .lLabbe et al.l I200E 



IWebb et all 120061 : iPapovich etaTT 120061: iReddy et all 
| 2006tl and very d eep spectra (jvan Dokkum et all 120041 : 
iKriek et"aT1l2006l ) has enabled the passively evolving sub- 
set of DRGs to be separated from the dusty subset, pro- 
viding more robust constraints on the physical proper- 
ties and evolution of massive > 10 11 M Q galaxies at 
z > 2. These studies have allowed much more accurate 
constraints to be placed on the epoch of for mation of the 
popul ation of massive galaxies at z ~ 2.3. IKriek et al.l 
(2006|) find that the near-infrared spectra of a small sam- 
ple of z > 2.3 DRGs are best fit by ages of < 1.3 — 1.4 
Gyr, suggesting that most of these systems started form- 
ing their stars between z ~ 3 and 4. 

While the formation epoch of the massive galaxies at 
high redshift is becoming increasingly better defined, it 
is still unclear what are the precurs ors of these mas- 
sive systems. Sub-mill i meter galaxie s (ISmail et al.lll998l: 
[ Chapman et~aTI 120031 : ICapak et all l2008t IDaddi et all 
l2008f) are convincing candidates since their intense star 
formation rates allow a substantial amount of mass to be 



assembled in a short period of time (see ICimatti et alJ 
12001 iMcCarthv et all 12004 L But sub -mm sources may 
not provide a complete sample of the precursors to 
high redshift massi ve and passive systems. Indeed, 
Shapl ev" et al.l (|2005f ) have shown that for reasonable as- 
sumptions the most massive objects in z ~ 2 rest-UV 
selected samples are likely to evolve into quiescent sys- 
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terns by z ~ 1.7, with space densities that are virtually 
identical to those o f the p assive population at that red- 
shift (|Cimatti et alj|2004f> . Objects may of course spend 
portions of their assembly history in both a SMG and 
LBG pha se. Our goal i n this section is to extend the 
efforts of Sh aplev et al.l (|2005h to quantify the fraction 
of massive galaxies at z ~ 2 — 3 that form a significant 
portion of their mass in an LBG phase at z > 3.5 and 
may be missed by current sub- mm galaxy surveys. 

We thus set forth to compute the number densities of 
massive objects in our B, V, and i'-drop samples. To 
derive the stellar mass functions, we group the dropouts 
by stellar mass in bins of AlogioM=1.0, spanning be- 
tween 10 7 5 and 10 115 M Q and then compute the num- 
ber density of dropouts in each mass bin, accounting for 
the variation of the effective volume with apparent mag- 
nitude. We also apply a correction to account for the 
galaxies which are not included in the SED-fitting anal- 
ysis due to being blended with foreground sources in the 
IRAC images. As we discussed in §5.1 (see also Figure 1), 
the fraction of sources that are blended varies with i' 775 
or Zg 50 apparent magnitude. We account for this slight 
bias by applying a magnitude-dependent correction: af- 
ter computing the number density of galaxies in a given 
stellar mass and apparent magnitude bin, we correct this 
density by dividing by the fraction of high-redshift ob- 
jects in that magnitude bin which were included in the 
SED-fitting analysis. In order to provide a consistent 
comparison of the mass function over z = 4 — 6, we im- 
pose an absolute magnitude limit of M1500 = —20. This 
ensures that each dropout sample is sufficiently complete 
to make a reliable measurement (see Figure 9). 

The resulting mass functions are shown in Figure 11. 
Also shown on the plot (purple datapoints) are the in- 
ferred densities if we include sources detected at 24/im 
in the mass function analysis. As expected, the 24/im 
detections slightly increase the number densities of the 
B and V-drops in the most massive bins, but generally 
do not increase them by more than the Poisson uncer- 
tainties in that bin. As discussed in §4.2, the inferred 
stellar mass uncertainties for individual galaxies become 
considerable for fainter sources with masses < 10 95 M . 
Although the average properties of sources less massive 
than this limit are reliable, we nevertheless choose to fo- 
cus our analysis on the more robustly determined objects 
contained in bins 10 95 M Q and above. 

The number densities derived here agree reasonably 
well with pre v ious m easurements. We plot the recent 
iMcLure et all (|2008l ) estimate of the z ~ 5 and 6 LBG 
stellar mass functions over our curves in Figure 1 1 (after 
adjusting their masses to be consistent with a Salpeter 
IMF) . In spite of their approximate way of determining 
the stellar mass function (the McLure et al. mass func- 
tion is derived by scaling the luminosity function by the 
average stellar mass to light ratio), the two estimates 
agree to within a factor of 2 over the mass range probed, 
with our determination predicting a larger density of the 
most massive galaxies. The results also are in reasonable 
agreement with the massive end of the stellar mass func- 
tion determ ined from the cos mological SPH simulations 
discussed in iNagamine et al.l (|2008l ) . The observed mass 
function becomes shallower than the model predictions 
at < 10 10 M Q . This may be due, in part, to obser- 
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Fig. 11. — Stellar Mass Function of B-drops (top panel), V-drops 
(middle panel), and i'-drops (bottom panel) with M1500 < —20. 
The green swath and associated solid black data points corresponds 
to the mass function determined from a large sample of dropouts in 
GOODS-S and GOODS-N. As discussed in §5.1, stellar masses are 
determined for individual galaxies by fitting the data against CB07 
population synthesis models. This plot is constructed from models 
with exponentially-decaying star formation histories with t=100 
Myr, solar metallicity, and a Salpeter IMF. The purple datapoints 
(slightly offset from the central mass bin for clarity) correspond to 
the densities inferred if 24/im-detected sources are included. The 
vertical dotted lines represent the mass limit above which the in- 
ferred masses are reliable for individual galaxies. The black dot- 
ted lines overlaid on the V- and i'-dropout mass functions (upper 
right and bo ttom panels) c o rrespo nd to the stellar mass function 
estimated in McLure et al. (2008) by multiplying the luminosity 
function by the average M/L ratio of the dropout population. The 
red dashed line overlaid on the i'-drop sample (bottom panel) cor- 
responds to the predictions of the z ~ 6 stellar mass function from 
the cosmological SPH simulations presented in Nagamine et al. 
(2008). The black solid curves overlaid in the upper right and 
bottom panel show the mass function determined in the B-drop 
sample (upper left panel). 
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vational incompleteness at low mass limits. However a 
similar "knee" in t he galaxy mass function is seen in the 
iDrorv et al.1 (|2005t ) data suggesting that the simulations 
m ay be underpredicting l ow-mass galaxies, as suggested 
bv lNagamine et all (|2008l ). 

In contrast to the weakly evolving mass function at z < 
1, the increase in stellar mass between z ~ 6 and 4 occurs 
over the entire stellar mass range considered. Between 
z ~ 6 and z ~ 4, we find a factor of 5-6 increase in 
density for the bins ranging between 10 9 5 and 10 10 5 Mq. 
There is some evidence that the number density grows 
more rapidly for the most massive galaxies. While the 
uncertainties are significant, the mass functions suggest 
that the number density in the bin centered at 10 11 M Q 
increases by a factor of 17 between z — 6 and z = 4, 
nearly three times greater than the growth in less massive 
bins. It seems likely that we are witnessing a formative 
period of massive galaxy formation. 

We finally turn to whether the observed number den- 
sities of massive LBGs at z > 2 can account for a sig- 
nificant component of the largely quiescent population 
of DRGs at z ~ 2. At z > 4, we find that the number 
density of dropouts with stellar mass greater than 10 11 
M© (excluding MlPS-detected sources) is ~ 2 x 10~ 5 
Mpc -3 (Figure 12). We find no overlap between these 
sources and the massive, sub-m m detected object s re- 
cently discovered in GOODS-N (|Daddi et al.ll200i . If 
these massive LBGs become quiescent on a timescale of 
less than 1 Gyr, then we would expect them to appear 
as red galaxies with little star formation at z ~ 2. In- 
deed if we evolve the colors and rest-UV luminosities of 
these massive LBGs forward in time to z ~ 2 assum- 
ing the simple exponential decay models considered in 
§6, we derive colors and optical apparent magnitudes 
that are consistent with t he observed properties of DRGs 
(jvan Dokkum et aT1l2006h . Comparing the observed den- 
sity of massive LBGs at z > 4 to the number density 
of quies cent 10 11 Mq sys tems at z ~ 2 — 3 (1.1 xlO -4 
Mpc~ 3 . lKrTek et al.ll2008h indicates that at least ~ 20% 
of the z ~ 2 — 3 population of quiescent 10 11 M galaxies 
formed a large fraction of their stellar mass at z > 3.5. 

Given the expectation that dust-enshrouded "bursts" 
of star formation seen in the sub-mm are largely re- 
sponsible for assembling massive galaxies b y z ~ 2 — 3 
(jCimatti et al.1 12004 iMcCarthv etaT1l2004h . it is inter- 
esting to extend our estimates of the number density 
of massive LBGs down to z ~ 2 in order to compute 
the total fraction of massive galaxies that were formed 
through objects pre sent in rest-UV selected samples. 
Shapl ev" et al.l (|2005f ) estimate that the number density 
of UV luminous galaxies with masses in excess of 10 11 
Mq at z ~ 2 is ~10~ 4 Mpc -3 . We make our own 
estimate by constructing a sample of rest-UV selected 
( 2 850 ^ 26.5) galaxies in GOODS-S with photometric 
redshfits between z — 2 and 2.4 using redshifts and 
photome try from the GOODS -MUSIC z 850 -band selected 
catalog (jGrazian et~aT] [20061 ). We follow the same pro- 
cedures described earlier for computing stellar masses, 
completeness as a function of magnitude, and effective 
volumes. The derived number densities are corrected for 
incompleteness (the sample is roughly 90% complete at 
Zg 50 = 26), resulting in final number density estimates of 
~ 9 x 10~ 5 Mpc -3 , very close to the measurement from 
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Fig. 12. — Evolving number density of massive galaxies at high 
redshift. Top: Solid black circles (and associated blue swath) show 
the evolution in the number density of actively forming galaxies 
more massive than 10 10 Mq derived in this paper. The green 
asterisk correspond to the total number density of galaxies more 
massiv e than 10 10 Mp) at z = 2.5 identified in a K s -band selected 
survey (Drory et al. 2005). The vertical dotted line demarcates the 
redshift range covered in this paper. Bottom: Same as top panel 
but for galaxies more massive than 10 Mq. The black circles 
give the number density of rest-UV selected galaxies more massive 
than 10 Mq; the upper datapoint at z = 4 corresponds to the 
number density inferred if the B-drops with 24/^m detections are 
truly at high redshift. The red asterisk shows the number density of 
galaxies more mas sive than 10 11 Mq that are quiescent at Z = 2.3 
(Kriek et al.ll2008l) . 



Shapl ev" 3t al.l (|2005l ). These densities suggest that at 
z ~ 2, the number density of UV- bright massive galax- 
ies is roughly comparable to the density of systems that 
have already become quiescent. 

Computing the contribution of rest-UV selected z > 2 
galaxies to the density of quiescent z ~ 2 systems is non- 
trivial as it depends strongly on the timescale over which 
star formation is quenched in the LBGs. The shorter the 
timescale, the larger the implied number density of qui- 
escent objects at z ~ 2. If we assume that massive LBGs 
join the red sequence over a period of 500 Myr, then we 
can estimate a quiescent number density by integrating 
over the growth in UV luminous massive galaxies over 
2 < z < 6 (Figure 12). Following this approach, we find 
that rest-UV selected galaxies would produce a quiescent 
galaxy density of 7xl0 -5 Mpc -3 by z ~ 2.5, increasing 
to 2xl0 -4 Mpc -3 at z ~ 2. If the quenching timescale 
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is longer (i.e., 1 Gyr), then the implied quiescent number 
density decreases to ~ 5 x 10" 5 Mpc" 3 at z ~ 2. Hence, 
these simple assumptions suggest that star formation in 
rest-UV selected z > 2 galaxies accounts for the assem- 
bly of at least 50% of the total quiescent population of 
massive galaxies at z ~ 2. 

It thus appears reasonable that a significant fraction 
of the passively-evolving subset of z ~ 2 — 3 galaxies 
with 10 11 Mg assembled their mass in UV luminous star 
formation episodes at higher redshifts. While it is pos- 
sible that these rest-UV selected galaxies would be de- 
tectable in deeper sub-mm observations, or alternatively 
that these galaxies may have undergone short-duration 
sub-mm luminous bursts of star formation in their past, 
these results appear to indicate that the ultr a-luminous 
sub-m m gal axies such as those presented in lCapak et al.1 
(2008) and iDaddi et ail (|2008f ) do not necessarily pro- 
vide the only route toward assembling massive, quiescent 
galaxies by z ~ 2. Constructing a complete sample of 
the precursors of this massive population of red galaxies 
must include both rest-UV and sub-mm samples. 

8. STELLAR MASS DENSITIES AT Z ~ 4, 5, AND6 

We now estimate the stellar mass densities implied by 
our dropout samples. By integrating the mass functions 
presented in Figure 11, we compute the stellar mass 
density brightward of the adopted rest-UV magnitude 
limit (taken to be Mi5oo=— 20 for each of the dropout 
samples). This absolute magnitude limit corresponds 
to apparent magnitude limits of i' 775 ~ 26.0 for the hy- 
drops, zgso ~ 26.4 for the V-drops, and zgso ~ 26.7 
for the i'-drops. Using the r=100(300) Myr exponen- 
tial decay models, we obtain 1.1(1.4) xlO 7 M Mpc -3 , 
3.7(4.9) xlO 6 M Q Mpc" 3 , and 1.6(2.3) xlO 6 M Mpc" 3 
for the B, V, and i'-drops, respectively (Figure 13). In- 
cluding all MlPS-detected dropouts sources increases the 
derived mass densities to 1.6(1.9) xlO 7 M Mpc" 3 for 
the B-drops and 4.7(6.2) xlO 6 M Q Mpc -3 for the V- 
drops with the r=100(300) Myr star formation history. 
The range of mass densities quoted above are consis- 
tent with o ur previous measu rements at z ~ 5 (5xl0 6 
M Q Mp c" 3 . iStark et allpOOl and z ~ 6 (2.5 xlO 6 M 
Mpc" 3 , lEvles et al.ll2007l) . 

As the total stellar mass density of galaxies at high 
redshift provides a unique constraint on the past star 
formation history of the universe, it is of interest to es- 
timate the contribution from galaxies that are fainter 
than the adopted rest-UV absolute magnitude limit of 
Mi5oo=-20. For the B and V-drop samples, this mag- 
nitude limit is brighter than allowed by the S/N limits 
of the data (see vertical lines in Figure 9). If we com- 
pute the mass density of the B-drop sample including all 
sources brighter than the 5a limit of M1500 — —19.1, we 
derive a mass density of 1.7(2.1) xlO 7 M Mpc -3 . Like- 
wise, integrating the stellar mass function of the V-drop 
sample including all sources brighter than the 5cr limit of 
Mi 500 ~ -19.5, we find a mass density of 4.4(5.9) xlO 6 
M Mpc" 3 . Both estimates exclude the contribution of 
MlPS-detected sources. 

In order to constrain the mass density of the popula- 
tion even fainter than these limits, we need to know the 
average stellar mass to UV luminosity ratio for galaxies 
below our detection thresholds. We can estimate this by 
extrapolating the trends present in the M*-Mi5oo relation 
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Fig. 13. — Stellar mass density of B, V, and i'-drops. The red 
rectangles demarcate the stellar mass densities derived from our 
dropout sample with an adopted magnitude limit of M1500 = —20. 
The circles (triangles) denote mass densities and associated Pois- 
son error determined using CB07 models with an exponentially- 
declining star formation history with a decay factor of r=100(300) 
Myr. See Figure 11 caption and §5.1 for further discussion of the 
assumptions used in deriving the stellar masses. Purple squares de- 
marcate rough estimates of the total stellar mass density of galax- 
ies brighter than M1500 = —10. These values are determined by 
extrapolating the average relationship between stellar mass and 
M1500 (Figure 9) to luminosities below the detection limits of the 
GOODS dataset. 

of Figure 9. Fitting a linear model to the observed trend 
for our B-drop sample, we find the following relation: 
logioM* = -0.26 - 0.47(Mi5oo). By extrapolating this 
function to fainter magnitudes and integrating the ob- 
served UV luminosity function (Bouwcns et al. 2 0071) to 



determine the number density of galaxies at each M1500, 
we can compute a rough estimate the contribution of 
typical faint UV sources to the total stellar mass density. 

Following this approach, we find that the entire pop- 
ulation of typical faint B-drops (e.g., those with average 
stellar masses in the magnitude range — 19 < Mi 50 o < 
—10) comprises ~ 49% of the total stellar mass density. 
This increases the total stellar mass density in B-drops 
to ~ 3.3(4.1) xlO 7 M Mpc" 3 for the t=100(300) Myr 
exponential decay models. Applying the same proce- 
dure to the V-drops (assuming the same relationship be- 
tween stellar mass and M1500 holds at z ~ 5), we derive 
total stellar mass densities of 1.0(1.4) xlO 7 M Mpc" 3 
for sources brighter than Mi5oo=-10. We refrain from 
making a similar extrapolation for the z'-drops given the 
larger uncertainties in this population's M*-Mi5oo rela- 
tion. 

9. SUMMARY 

We have studied the growth of galaxies over 4 < z < 6 
with the goal of improving our understanding of galaxy 
evolution in the first 2 billion years of cosmic history. 
Using the deep HST data available in the two GOODS 
fields, we compiled samples of 2819 B-dropouts, 615 V- 
dropouts, and 166 i'-dropouts. The superior quality of 
the multiwavelength data in GOODS (including opti- 
cal data from HST spanning four filters, near-infrared 
data from the VLT and Subaru, and mid-infrared data 
from Spitzer) greatly helps identify foreground galaxies 
and stars which lurk in Lyman break dropout samples. 
Great care was taken to identify and remove these ob- 
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jects from our sample, using HST morphologies to excise 
point sources likely to be bright stellar contaminants and 
broadband SEDs spanning from the optical to the in- 
frared to remove low-redshift galaxies. 

Our final high-redshift sample consists of 2443 B- 
dropouts, 506 V-dropouts, and 137 i'-dropouts. From 
this sample of galaxies, we constructed sub-samples of 
dropouts which are sufficiently unblended with nearby 
foreground objects in the deep Spitzer data to allow accu- 
rate photometry. In general, this corresponds to > 35% 
of the total samples. We used the CB07 population syn- 
thesis models to estimate stellar masses and ages for this 
"Spitzer clean" subset. We compared the inferred prop- 
erties to those derived from BC03 models, finding general 
agreement for solar metallicity models when the redshift 
is above 5. Below this redshift, the mid-IR data begin to 
probe the rest-frame near-infrared (which is affected by 
TP-AGB stars), and the CB07 models therefore return 
masses which are typically 10% lower than those from 
BC03 models. 

Using this large database of stellar masses and rest- 
UV properties of B, V, and i'-dropouts, we examine the 
evolution of star-forming galaxies between z ~ 4 and 6. 
We summarize our main results below. 

1. We find that the typical stellar masses and ages of 
galaxies of a fixed UV luminosity (uncorrected for dust 
extinction) do not evolve strongly between z ~ 6, 5, and 
4. This argues against the notion that the B-drops are 
made up of galaxies that have been steadily growing at 
fixed UV luminosity since z ~ 6. Instead, these data 
suggest that each successive epoch studied in this survey 
(e.g., z ~ 4, 5, and 6) is likely dominated by galaxies that 
have only recently (within the last ~ 300 Myr) emerged 
at their present luminosity. This is consistent with the 
observed dec line in the UV luminosity function between 
z ~ 4 and 6 (jBouwens et al.ll2007h . 

2. The evolution in the M^-Misoo relation and the UV 
luminosity function is potentially difficult to reconcile if 
the galaxies observed at z > 4 assembled their stars with 
a constant star formation rate. In this scenario, the in- 
ferred ages of galaxies at one epoch are sufficiently large 
to violate the required decline in the UV LF and the lack 
of strong evolution in the M^-Misoo- However if lower 
metallicity templates are adopted, the ages are lowered 
and the discrepancy reduced. Deeper infrared photome- 
try will be of significant benefit in improving constraints 
on the inferred stellar populations. 

We also considered a steady growth scenario in which 
star formation is becoming steadily more vigorous with 
increasing cosmic time betwe en z ~ 6 and z ~ 4 , as 
predicted in the simulations of iFinlator et al.l (|2006f ). In 
this case, the inferred ages of the galaxies remain large, 
but since their precursors are fainter in the past, they 
don't violate the observed evolution in the UV LF and 
M*-Mi5oo relation. The only problem that potentially 
stands in the way of this model is that it is inconsistent 
with the short duty cycles imp lied by measure ments of 
the clustering of z > 4 LBGs (|Lee et al.l [2001) . Addi- 
tional spectroscopy and new multiwavelength imaging in 
independent fields are required to confirm these cluster- 
ing measurements. 

3. We argue that episodic star formation may provide 
the best fit to the wide range of observations considered 



in this paper. In this picture, each epoch is dominated 
by recently emerged sources, which remain luminous for 
no more than ~ 500 Myr on average. If the UV lu- 
minous periods become more intense with cosmic time 
between z ~ 6 and 4, then this scenario could explain 
the evolution of the M+-M1500 relation and the UV LF. 
If this is the dominant mode of star formation at z > 4, 
then there should be a significant population of relatively 
massive, quiescent sources that are in between star for- 
mation episodes. Future efforts at z ~ 4 should help 
clarify whether this is the case. 

4. The stellar mass function of UV luminous sources 
grows significantly from z ss 6 to z ss 4, typically in- 
creasing by more than a factor of six in number den- 
sity over the entire stellar mass range considered. There 
is some evidence that the growth in the most massive 
LBGs (M* ~ 10 11±0 5 ) is a factor of three greater than 
in less massive systems, implying that this redshift range 
likely corresponds to the beginning of massive (10 11 M ) 
galaxy formation era. 

5. The rapid growth in the assembly of massive galax- 
ies at z > 4 indicates that at least ~ 20% of the z ~ 2 — 3 
quiescent massive galaxies assembled their stellar mass 
at z > 3.5. Fitting the increase in number density of 
UV luminous, massive galaxies between z ~ 6 and 2, 
we estimate that it is feasible that more than 50% of 
quiescent z ~ 2 — 3 massive galaxies were assembled in 
higher redshift LBGs. These results imply that the high 
star formation rate, dust-obscured sub-mm galaxies are 
not the only route toward assembling quiescent 10 11 M Q 
systems by z ~ 2 — 3. 

6. Using the CB07 models with solar metallicity and 
a exponentially declining star formation history with a 
decay factor of r=100(300) Myr, we derive the stellar 
mass density of rest-UV selected dropouts brighter than 
M1500 = -20: p± ~ 1.1(1.4) xlO 7 M© Mpc" 3 at z ~ 4, 
3.7(4.9) xlO 6 M Mpc" 3 at z ~ 5, and 1.6(2.3) xlO 6 
Mq Mpc -3 at z ~ 6. These measurements are certain 
to underpredict the total stellar mass density since they 
do not include the contribution from sources that either 
do not satisfy are color selection criteria or lie faintward 
of the absolute magnitude limit. By extrapolating the 
M*-Mi5oo relation from Figure 9 to M 1R n n=-10 and in- 
tegrat ing the UV luminosity functions of iBouwens et al.1 
(|2007h . we estimate that the population of faint galaxies 
may increase the mass density estimates quoted above 
by a factor of 3-5 x . 
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